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Abstract 



We directly fit the experimentally measured energy dependence of the average value 
of 1 — Thrust, <1 — T>, over the e + e~ centre-of-mass energy range Q = 14 — 172 GeV 
to the QCD expectation obtained by integrating up the evolution equation for 
d<l — T>/d\og Q in terms of <1 — T>. We fit for A^L, uncalculated C(a|) and higher 
perturbative corrections parameterized by the scheme invariant p2, and the parameter 
A which characterizes the magnitude of the leading 1/Q power corrections anticipated 
for <1 - T>. A 3-parameter fit yields = 245±?° MeV, p 2 = -16 =F H and 

A = 0.27i°;io GeV, equivalent to a s (M z ) = 0.1194 ± 0.0014. In this approach, there 
is no error associated with the choice of the renormalization scale u. 



For several e + e QCD observables we now have experimental measurements Jl|, |2| span- 
ning the centre-of-mass energy range from the lowest PETRA energy Q = 14 GeV through 
LEP-I at Q = 91 GeV up to LEP-II at Q = 172 GeV. An ideal observable for testing QCD 
is the average value of 1 — Thrust =<1 — T>, where the Thrust is defined to be, 

rp Ek\Pk-n\ 

1 = max , (1) 



with the sum running over all particles in the event. The thrust axis n is varied to maximize 
the thrust. Thrust describes the jettiness of the event such that T = 1 for events with 
two back-to-back particles and T = 1/2 for completely spherical events. Since T is fully 
inclusive, the averaging means that it is free of the large kinematical logarithms which afflict 
distributions in jet observables close to the two-jet region. 

If we consider the observable R(Q) =<1 — T>/1.05 then we have a perturbation series 
and leading power correction of the form, 

R(Q) =a + r x a 2 + r 2 a 3 + . . . + — (1 + \ x a + X 2 a 2 + ...), (2) 

where a = as^/n denotes the renormalization group improved coupling. Note that the 
normalization is simply such that the perturbative expansion begins with unit coefficient. 



In the MS scheme with \i = Q and Nf = 5 active quark flavours the next-to-leading order 
(NLO) coefficient is T\ = 9.70 H]. The next-to-next-to-leading order (NNLO) coefficient 
T2 is unknown. Other uncalculated higher order corrections and genuine non-perturbative 
corrections are included by the phenomenological 1/Q power corrections 0, |, || . In , the 
power corrections up to A 2 are calculated in terms of a scale /// representing the transition 

1 



between the perturbative and non-perturbative regimes. To extract a s {Mz) from the data, 
we just truncate the perturbative series for a given renormalization scale \i = xQ (which 
is typically x — 1), and a given value of pi (typically 2 GeV). Then, by comparing with 
experimental data, we solve for a. A recent analysis for <1 — T> using this approach 
finds, 

a s {M z ) = 0.1204 ± 0.0013 Iffi j£S (3) 

(with a xVd.o.f of 42.6/24) where the first error is purely experimental. The second and 
third errors come from varying the theoretical input parameters, first allowing x between 
0.5 and 2 and second for \ii in the range 1 — 3 GeV. Clearly the estimate of the theoretical 
error is dominated by the renormalization scale uncertainty. 

Alternatively, we may avoid the renormalization scale entirely and, by differentiating 
Eq. H] with respect to Q and using the renormalization group equation for the running of a, 
directly write an expression for the running of R{Q) with Q in terms of R{Q) itself ]7], || |, [TI| , 

— — = -bR 2 {l + cR + p 2 R 2 + ...) + K Q R- c l b e- 1 ' hR {l + K 1 R+...) + ... 
dlogQ 

= -bp(R). (4) 
Here b and c are the first two universal terms of the QCD beta-function, 

t - (5) 



153 - 19JV) ... 

c = -. (6) 

126 v ; 

The quantity, 

p 2 = r 2 + c 2 - nc - rl, (7) 



is a renormalization scheme and renormalization scale (RS) invariant combination of the 



NLO and NNLO perturbation series and beta-function coefficients with, in the MS scheme, 

77139 - 15099iV f + 325N 2 

C2 = msT (8) 

Since the NNLO r 2 is unknown, so is pi- As we will see later, the coefficient Kq is directly 
related to the coefficient A of the 1/Q power corrections in eq. (0). 

Since R{Q) and dR/d\ogQ are both observables one could in principle directly fit eq. (f|) 
to the data and thus constrain the unknown coefficients p 2 and K Q . At asymptotic energies 
all observables run universally according to, 

-^5- = -&#(! + eft), (9) 
d log Q 

and one could see how close the data are to this evolution equation. Given the error bars 
of the data and the separation in Q of the different experiments it is preferable, however, 
to integrate up eq. (f|) using asymptotic freedom (R(Q) — > as Q — ► oo) as a boundary 
condition. In this way one obtains 0, 

i + clogf-^-)=Mogf^-)- f R dx (—}-+ 1 ), (10) 

R \l + cRJ \ArJ Jo V P( x ) x 2 (l + cx)J v ; 

where A# is a constant of integration. By comparing with the Q — ► oo behaviour of eq. (0) 
one can deduce that |J, 

A« = e^( y ) Ajrs, (11) 
where r = rf ts (p = Q). Evaluating this for < 1 — T > yields A# = 14.4A^r. 

If we assume that the right-hand side of eq. (f§) is adequately parameterized by, 

- bp(R) = -bR 2 {l + cR + p 2 R 2 ) + K R~ c/b e- 1/bR , (12) 



we can then insert this form into eq. (|10"D and by (numerically) solving the transcendental 
equation, perform fits of p 2 , K and to the data R(Q). The fitted K can then be 
converted into an estimate of the parameter A in eq. (|2|) by differentiating eq. (0) with 
respect to log Q and comparing terms. We find, 

fh\ c,b 

\ = -K e"> l^-j A m . (13) 

In Fig. |I] we show the fit to the data obtained by setting p 2 = A = 0. This corresponds to 
the universal running of the observable given in eq. (g). The fitted value is A§L = 266 MeV 
with a % 2 /d.o.f = 81.7/32. The dashed lines show the effect of varying Aj^g by ±30 MeV 
around the fitted value. We clearly see that the data is falling much too quickly with 
increasing Q for the asymptotic behaviour to have set in at these scales. The data favours a 
more steeply falling evolution which could be caused by either higher order corrections with 
a positive p 2 , or power corrections with non-zero Kq. We therefore perform a 3-parameter 
fit allowing p 2 , K and A-^ to vary independently which is shown in Fig. ^. The minimum 
X 2 fit is acceptable (x 2 /d.o.f = 40.4/30) and estimating an error by allowing \ 2 within 1 of 
the minimum gives, 

Ag_ = 245^ 2 ? MeV with p 2 = -16 T H and A = 0.27j£Jjj GeV. 

These values of p 2 and A are reasonably small, thereby lending support to our critical as- 
sumption that the evolution equation could be parameterized in this way. We can convert 
the extracted value of Ajj^ into as(Mz) using, 



a s (M z 



Mog(QVAfe) 



i 2clog(log(Q 2 /A^)) 



&log(Q7A§^ 



2 / A 2 _ 

(14) 



and find, 



a s (M z ) = 0.1194 ±0.0014. 



(15) 



The fitted p 2 value may also be converted, if desired, into an estimate of r 2 (/i — Q) — 
89 ± 11. Note that this large value of r 2 is almost entirely due to the renormalization group 
predictable r\ + cr\ piece relating p 2 and r 2 (see eq. (|7|)). 

We see that our central value is remarkably close to that obtained by [[|. The main 
difference is in how the errors are determined. In our approach, the errors are estimated by 
allowing the uncalculated higher orders to be fitted by the data and the data prefers these 
to be small. As higher order corrections become known, the new RS-invariant terms, p 2 , P3 
etc., can be incorporated and the fit refined. 

With such an accurate value of as(Mz), we should expect that applying this approach 
to other variables should yield consistent results. Unfortunately, the method described here 
relies on having reliable data over a reasonable range of Q values. Another variable for 
which a wide range of data has been accumulated is the heavy jet mass, which is obtained 
by assigning the particles to one of two hemispheres Hi and if 2 according to whether pk ■ n 
is positive or negative, and finding the maximum scaled invariant mass of the hemispheres, 



Here R(Q) =<m 2 H /s > /l. 05 while r\ = 4.52 If we repeat the same analysis for this 
variable, we find that a one parameter fit with p 2 = A = gives a very poor fit, x 2 /d.o.f = 
213/29 and = 368 MeV. As seen in Fig. ^ the data evolves much faster than the QCD 




(16) 
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prediction. However, allowing both p 2 and A to vary while using the value of Ajjg = 245 MeV 
obtained from the <1 — T> analysis gives a much more satisfactory description.^] Here, 
X 2 /d.o.f = 40.3/27 while p 2 = 13 and A = 0.11 GeV are sufficiently small to support our 
choice of parameterization. 

It is clear that the only certain way to reduce theoretical errors is to compute the NNLO 
coefficients p2- However, fitting directly to the energy dependence of the data as we have done 
enables us to estimate them together with possible power corrections. The renormalization 
group scale-dependent logarithms are automatically resummed to all orders on integrating 
eq. (H) and do not add a spurious extra large uncertainty in the extraction of Aj^g (or 
equivalently a s (M z )). As a result, we have obtained a value of the strong coupling constant 



that is both competititive with and in agreement with the current world average |TT[ of a 
wide range of data, 

a s (M z ) = 0.119 ±0.002. 
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1 Unfortunately, in a three parameter fit, p2 and A trade off against each other and drive pi and A to 
unacceptably large values where we would have no reason to believe that the parameterization is adequate. 
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Figure 1: The average value of 1 — T obtained experimentally [JTJ, compared with the 
QCD expectation of eq. fllPp. The solid line shows the fit to the data with p 2 = A = 
corresponding to = 266 MeV. The long-dashed and short-dashed lines show the effect 
of altering A^ by ±30 MeV. 
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Figure 2: The average value of 1 — T obtained experimentally [||, compared with the QCD 
expectation of eq. (0). The dashed line shows the fit to the data with p 2 = A = while the 
result of the three parameter fit (to Ajjg, p 2 and A) is shown as a solid line. 
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Figure 3: The average value of the heavy jet mass obtained experimentally |], |2[ compared 
with the QCD expectation of eq. flTOl). The dashed line shows the fit to the data with 
P2 = X = while the result of the two parameter fit (to p2 and A) using the value of 
obtained from the three parameter fit to <1 — T> is shown as a solid line. 
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